Liquid crystalline states for two-dimensional electrons in strong magnetic fields 
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Based on the Kosterlitz-Thouless-Halperin-Nelson- Young (KTHNY) theory of two-dimensional 
melting and the analogy between Laughlin states and the two-dimensional one-component plasma 
(2DOCP), we investigate the possibility of liquid crystalline states in a single Landau level (LL). 
We introduce many-body trial wavefunctions that are translationally invariant but posess 2-fold (i.e. 
nematic ), 4-fold (tetratic) or 6-fold (hexatic) broken rotational symmetry at respective filling factors 
v = 1/3, 1/5 and 1/7 of the valence LL. We find that the above liquid crystalline states exhibit 
a soft charge density wave (CDW) which underlies the translationally invariant state but which 
is destroyed by quantum fiuctuations. By means of Monte Carlo (MC) simulations, we determine 
that, for a considerable variety of interaction potentials, the anisotropic states are energetically 
unfavorable for the lowest and first excited LL's (with index L = 0, 1), whereas the nematic is 
favorable at the second excited LL (L — 2). 
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I. INTRODUCTION 

In 1983 Laughlin [1] introduced his famous trial wave 
function 
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to describe the fractional quantum Hall efFect(FQHE) 
states [2-5] for filling factors v = 1/m oi the lowest Lan- 
dau level (LLL), where m is an odd integer. Immediately 
after this discovery, many attempts were done to compare 
the stability of these states against other known ground 
states, typically Wigner crystal (WC) states [6-8]. At ab- 
solute zero (T — 0) the current theoretical understanding 
is that WC states are favorable for filling factors smaller 
than a critical value i^c — 1/6.5 [7,8]. For larger filling 
factors of the LLL, the electrons are believed to form a 
quantum liquid state with Laughlin wave function being 
an excellent choice for v = l/m (with m = 1,3,5) [9]. 
Because of its translational and rotational invariance, 
Laughlin's wave function can be used to describe a liquid 
state of the electrons in the LLL, as can be seen by writ- 
ing as a classical distribution function [10,11]: 
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and V is the potential energy of a classical two- 
dimensional one-component-plasma (2D0CP) system. 



Using the formal analogy between the Laughlin wave 
function and the 2D0CP we can identify a dimension- 
less coupling constant, F = /3e^ = /{ksT) = 2m. An 
equlibrium state of the 2D0CP is entirely characterized 
by F and the freezing transition in this case was located 
at F « 140 [10]. Employing the analogy between the 
temperature of the classical plasma and the filling fac- 
tor of the LLL, we should expect a freezing transition as 
we decrease the electronic filling factor in the quantum 
Hall regime. Because of the different quantum nature 
of the electronic correlations in the FQHE, it was found 
that such a system is a Laughlin liquid for filling factors 
V = 1/3 and 1/5, but becomes crystal for filling factors 
smaller than « 1/6.5 (this value is about an order of 
magnitude larger than that deduced from the classical 
2D0CP analogy). 

It is feasible that, in analogy to the classical freez- 
ing transition realized by cooling down a 2D0CP, the 
transition to a solid (WC) state obtained by reducing 
the filling factor in the electron case may be interpreted 
as a topological Kosterlitz-Thouless-type transition [12]. 
This would be the correlated electron system counter- 
part of the wellknown 2D melting problem. Although 
the 2D melting is not fully understood, an elegant and 
reliable theory of melting has been proposed in the 1970's 
by Kosterlitz, Thouless, Halperin, Nelson and Young 
(KTHNY) [12-14]. The KTHNY theory predicts that 
an intermediate third phase called hexatic, will exist be- 
tween the hexagonal solid and the liquid phases in a cer- 
tain portion of the phase diagram (perhaps in a some- 
what narrow range of temperatures). In the liquid phase 
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there is no long-range translational or rotational order 
(the system is both translationally and rotationally in- 
variant). In the solid phase the system has quasi- long- 
range translational and true long-range rotational order. 
The hexatic phase in the KTHNY theory is thought to 
have no true long-range translational order, but does re- 
tain quasi-long-range oricntational order (the system is 
translationally invariant, but not rotationally invariant 
at least for short distances). The intermediate hexatic 
phase is often considered most important since it has a 
symmetry intermediate between the hexagonal solid and 
the liquid. 

Recent experiments in very high mobility {n ^ 10^ 
m/Vs) GaAs/Ala:Gai_a;As heterostructures have shown 
a variety of low temperature phases with exotic proper- 
tics. Since 1999 it has been known that in transitional re- 
gions between QH plateaus for high LL's (with LL index 
L > 2) either a smectic or nematic phase exists [15-20]. 
In fact, one of us calculated to a reasonable accuracy the 
anisotropic-isotropic transition temperature as a topolog- 
ical process [17]. In 2002 a melting transition from the 
WC state to a FQHE-like state was observed at ca. 130 
mK [21] and speculation mounted to suggest that possi- 
bly this transition occurs to a hexatic mesophase [22]. 

On this grounds we investigate the possibility of vari- 
ous liquid crystalline mesophases in a partially filled LL. 
Given that two-dimensional liquid crystals may posses 
different forms of rotational group symmetry, wc select 
a set of possible candidates, having C2 (nematic), C4 
(tetratic), and Cq (hexatic) rotational group symmetry 
(note that in principle higher symmetry groups are also 
possible for a liquid crystal, e.g. a liquid quasicrystal 
with a Cio symmetry — we have not explored, however 
such possibilities in this paper [23]). Our results indi- 
cate that the states studied exhibit a soft charge density 
wave (CDW) which underlies the translationally invari- 
ant state but which is destroyed by quantiim fluctuations. 
We perform Monte Carlo (MC) simulations and deter- 
mine that, for a wide range of interactions the anisotropic 
states are energetically unfavorable for the lowest and 
first excited LL's (with index L — 0, 1), whereas the ne- 
matic is favorable at the second excited LL (L = 2). 

In Sec. II we describe the types of states that were 
considered for our calculations. Section III presents the 
types of interaction potential considered and explains the 
methods used to calculate the properties of the system. 
Section IV contains the results obtained and a discussion 
of their meaning. The underlying soft CDW is discussed 
in Sec. V. Finally, the conclusions are presented in Sec. 
VI. 



II. LIQUID CRYSTAL STATES 

In this paper we consider liquid crystalline phases with 
no translational order but with quasi-long-range orien- 



tational order with various rotational symmetry groups 
C2 , C4 , and Ce ; corresponding to a nematic, tetratic and 
hexatic phase respectively. There are some basic require- 
ments on how we construct these states: (i) the states 
must obey Fermi statistics, i.e. they must have odd par- 
ity under the exchange of any pair of electrons; (ii) the 
states must be translationally invariant (at least far away 
from the boundaries of the system in case of a finite num- 
ber of electrons); (Hi) there must be a broken rotational 
symmetry belonging to the proper symmetry group; (iv) 
the states must belong to a single LL to avoid the large 
cyclotron energy cost Uwc = heB/nie, where B is the 
magnetic field, and e and rrig are the electron charge and 
mass respectively (also note that as we will show later, 
various properties at any LL can be readily obtained from 
properties calculated in the LLL). 

A class of such wave functions satisfying all these re- 
quirements are the so-called broken-rotational-symmetry 
(BRS) wave fimctions [18-20, 22, 24] that are systemat- 
ically constructed by properly splitting the zeros of the 
Laughlin liquid state [in essence, the idea is to place the 
vortices that perform the composite fcrmion (CF) trans- 
formation [4,5] around the location of the electron, rather 
"on top" of them]. Let us consider the Laughlin wave- 
function as given in Eq. (1), where Z). = Xk + iyk is fc-th 
electron position in the xy-plane in complex notation, 
and lo = [h/{eB)]^/'^ is the magnetic length. This wave 
function represents a gaped, uniform and isotropic liquid, 
and is an excellent description of a liquid state at filling 
factor u = 1, 1/3 and 1/5 of the LLL (for u = 1/7, the 
WC state prevails, see previous discussion, and Ref. [22]). 

To build a liquid crystal (BRS) state out of the liq- 
uid states we split the zeros of the wave function in a 
way that conserves the anti-symmetry (Fermi statistics) 
and translational invariance, but breaks the rotatational 
invariance of the wave function. This is done by intro- 
ducing a prefered set of directions [18 20,22,24] into the 
wave function creating a degree of anisotropy. A gen- 
eralized liquid crystal wave function for a filling factor 
u = 1/m can then be easily written as: 
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where the complex directors are distributed in pairs 
of opposite value in the complex plane (to satisfy Fermi 
statistics). In this paper we focus on the states with the 
highest level of discrete symmetry possible at each filling 
factor, which is set by distributing the symmetrically 
in a circle around the origin: 

a^ = ae'2"(''-i)/('"-i), /x e {1, 2, . . . , (m - 1)} . (4) 
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Without loss of generality a can be taken to be real. 
The wavefunction in Eq. (3) represents a homogeneous 
liquid crystalline state at filling factors v = 1/m, is anti- 
symmetric, lies entirely in the LLL, and is smoothly con- 
nected to the isotropic Laughlin state for a = 0. 



III. INTERACTION POTENTIALS AND MONTE 
CARLO SIMULATION 

For our simulations we consider N electrons in a charge 
neutralizing background. When considering the quan- 
tum Hamiltonian H = K + V, the strong magnetic field 
quantizes the kinetic energy K so that single-LL wave- 
functions have a constant (and thus irrelevant) kinetic 
energy, {K)/N. The only relevant contribution comes, 
therefore, from the total potential energy operator 



V = Vee+ Veb + Vbb , 



(5) 



consisting of electron-electron, electron-background and 
background-background interactions. 

It has been a common practice to work on the sur- 
face of a sphere [5] in order to minimize boundary ef- 
fects in the finite-size computations. However, due to 
the anisotropic nature of the states under consideration 
this scheme would produce significant problems due to 
the need to have topological defects at the "poles" of the 
sphere. We therefore work on a simpler disk geometry, 
where the neutralizing positive background has a uniform 
density po — v / {2t:11) and is spread over a disk of radius 
iJjY = Iq{2N/vY/'^ with an area ri^r = ttR\. 

Our goal is to thoroughly investigate the possibility 
of a liquid crystal state in the LLL for electrons inter- 
acting not only with the usual bare Coulomb potential 
vc{f 12) = e^/ (e'"i2) but also for a variety of other reason- 
able effective potentials that take into consideration the 
finite thickness of the quasi-2D electron layer. As pre- 
viously shown by Zhang and Das Sarma (ZDS) [27], the 
electron-electron interaction in a quasi-2D system can be 
written as: 
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where ri2 is the 2D distance separating the two electrons, 
e is the average background dielectric constant, Jq is the 
Bessel function of zeroth order, and 6 is a parameter re- 
lated to the finite thickness of the 2D layer (if we define 
the average thickness as Z, then b = 2>/ Z). In addition, 
we also consider two other interaction potentials: 
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The two model potentials include the thickness effect 
phenomenologically [27] through the length parameter, 
\ = Z /2 = l.b/b. All the above potentials have the same 
Coulomb behavior for large ri2, but differ from the bare 
Coulomb potential for small ri2- 

To consider the zero-temperature stability of the liq- 
uid crystal states of Eq. (3) with respect to the uniform 
isotropic liquid state counterparts, we performed exten- 
sive MC simulations in order to compute the energy and 
other quantities for the four different interaction poten- 
tials. Since the potentials involved are merely single- and 
two-body interactions, we need to accurately determine 
all single- and double-particle distribution functions, i.e. 

the density p(r) = ^X]i!=i ^i^i ^ ^)'^^ ^^'^ the pair corre- 
lation function g{r 12), respectively. The determination of 
such functions allows an accurate determination of all po- 
tential energies in the N ^ 00 thermodynamic limit [28] . 

By definition, the pair correlation function, g{ri2) is 
the conditional probability [normalized so that g{oo) = 1] 
to find an electron at position r" given that another elec- 
tron is found at position r' = r" — ri2: 
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where po = i^/{2ttIq) is the average bulk electron density. 
It is also useful to define the static structure factor /S'(q), 
which is given by the 2D Fourier transform of g{ri2)'- 
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Note that, because of the anisotropy of the wave function, 
both functions are explicitly angle-dependent: .g(ri2) = 
9{'i"i2,6) and S{q) = S{q,9q) for a ^ 0. It is also worth 
noting that the charge neutrality sum rule guarantees 
that S'(q) cx ^2 for g ^ [10,25,26]. 

In the thermodynamic limit, the ground state correla- 
tion energy per particle can be easily computed from [28]: 

= ^ = ? / ^'""^^ ''^''''^ t^^'"^'^ ~ ^] ' 

where v{ri2) can have any reasonable form, in particular 
it can take the form of any of the potentials shown in 
Eqs. (6,7). Because the interaction potentials are cen- 
trally symmetric, the above formula can be rewritten in 
the simpler form: 

POO 

Ea = ^(27r) / dri2 ri2 ^;(ri2) [g{n2) - 1] , (11) 
^ Jo 

where g{r 12) is the angle-averaged pair distribution func- 
tion: 
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For specific cases [as for the VzDsij'12) potential, which 
has strongly oscillatory behavior in real space, making 
the numerical calculations very unstable and precarious] 
a corresponding formula that uses the static structure 
factor was employed: 

11/"°° 

In this case v{q) is the 2D Fourier transform of the inter- 
action potential, and we also define the angle-avcragcd 
static structure factor S{q) = J^^ d0q/{2Tr) S{q). The 
use of the static structure factor has the added advan- 
tage of allowing the calcuation of the correlation energies 
in all LL's from a single determination of the pair cor- 
relation function in the LLL [18-20]: 
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where Ll{x) are Laguerre polynomials and L corre- 
sponds to the LL index. 

As in any MC calculation using the the Metropolis 
algorithm [29] , the expectation value of any [position de- 
pendent, e.g. p{t)\ operator can be computed by averag- 
ing the local value of the operator over a large number of 
electronic configurations generated from the probability 
distribution P oc [^'"y^l^. In a MC attempt, one eletron 
is moved to a new position Y^riai = + A^, where 
is a random vector in some domain. If the probability 
ratio, P{rtrial)/P{^i) is larger than a random number 
uniformly distributed in the [0,1] range then the move 
is accepted and we let let rj+i = r trial, otherwise the 
move is rejected and rj+i = r^. We adjust the size of 
the domain over which A^'s vary so that about half of 
the attempted moves are accepted. Following standard 
practice, we denote a MC step (MCS) a sequence of steps 
described above so that every electron in the system has 
attempted a move (and about half succeed). After a MCS 
the system is in a state essentially uncorrelated to the 
previous one and averages are computed for the desired 
operators [30] The results we report were obtained after 
discarding 100,000 "thermalization" MCS's and using be- 
tween 2 X 10^ and 4 x 10^ MCS's for averaging purposes 
on systems of 200-400 electrons. 



IV. MONTE CARLO RESULTS AND 
DISCUSSION 

By using MC methods we studied the possibility of a 
liquid crystal state in the LLL for the leading candidate 
states at filling factors, u = 1/3, 1/5 and 1/7. A trial 
wave function as in Eq. (3) was considered and various 



properties were analyzed as function of the anisotropic 
parameter a. Various interaction potentials were consid- 
ered for the computation of the correlation energies [see 
Eqs. (6, 7, 10, 11, 13, 14)], all have in common the fact 
that they incorporate the effects of finite layer thickness 
into the quasi-2D electronic system and are essentially 
identical to Coulomb's for large distances. This choice 
is motivated by the wellknown fact that the finite layer 
thickness of a real 2D system leads to a weakening and 
eventual collapse of the FQHE [31]. Therefore, when the 
finite layer thickness (parameter A) increases as to be- 
come larger than the magnetic length, the short-range 
part of the Coulomb interaction softens and as a result 
the isotropic FQHE liquid state may become unstable 
with respect to another state of different nature (a pos- 
sible new candidate can be the liquid crystal state con- 
sidered here, and/or a Wigner crystal). 
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FIG. 1. Angle-averaged single-particle density, p(r), for 

N = 196 electrons and filling factors i/ = 1/3, 1/5 and 1/7. 
We show the results for the isotropic cases (o = 0) and, for a 
large a (the oscillations observed in this case arc discussed in 
Sec. V). Here r is the distance from the center of the disk. 

In Fig. 1 we show a plot of the angle-averaged single- 
particle density p(r) for states oi N = 196 electrons and 
fining factors of 1/ = 1/3, 1/5 and 1/7. The existence, 
for a = of a large region around the center of the disk 
(r = 0) with constant density is an indication that there 
is bulk-like behavior [28]. Results for moderate values of 
a are similar to those for a = 0. For larger a an appar- 
ent density fluctuation propagates from the edges to the 
center making it very difficult to identify a "bulk" region. 
The existence of this density fluctuation is discussed in 
detail in Sec. V. We found that values of a acceptable 
for the purposes of calculating bulk-like properties in rea- 
sonably sized systems are as follows: a < 3 for i/ = 1/3, 
a < 4 for both v = 1/5 and 1/7 respectively. 

In order to compare the energy of the isotropic Laugh- 
lin liquid state with that of an anisotropic liquid crystal 
state, we first need an accurate computation of the pair 
distribution function in terms of the parameter a. For 
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the smallest a's, a number of iV = 196 electrons was 
sufficient to give a very accurate pair distribution func- 
tion, whereas as many as 400 electrons were used when 
as became large as to induce sizeable oscillations in the 
density. Figure 2 shows results for the pair distribution 



function, g{r), for the: v = 1/3, a = 2 nematic, v = 1/5, 
a = 3 tetratic, and v — 1/7, a = 3 hexatic. Each MC 
simulation involved 4 x 10^ MCS's and ca. 400 electrons. 
Figure 3 shows the corresponding static structure factors 
S'(q) obtained from g{r) using Eq. (9). 




FIG. 2. Pair correlation function g{r) for u — 1/3, a — 2 (left panel), v — 1/5, a — 3 (center panel), u = 1/7, a = 3 (right 
panel) . Note the discrete rotational symmetry of each state. 



v=l/a v=I/5 v=l/7 




FIG. 3. Static structure factor ^(q) for u — 1/3, a = 2 (left panel), u = 1/5, a = 3 (center panel), v = 1/7, a — 3 (right 
panel). Note the discrete rotational symmetry of each state. 



Since the angle-averaged 3(^12) is sufficient for the determination of the energy, we averaged it (at significant savings 
in computer time) for various combinations of filling factor v, and anisotropy parameter a. Figure 4 shows some of 
our results for 196-400 electrons. 
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FIG. 4. Angle- averaged pair correlation function (?(ri2) for 1/ = 1/3, = 1/5 and v = 1/7. 
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At all filling factors that wc considered, we noted that 
g{i^i2) changes very little when parameter a is small (e.g. 
1). Only for larger a's (> 2) sizeable changes take effect. 
In view of this behavior, wc anticipate that the energy 
differences between the isotropic liquid state (a = 0) and 
the anisotropic liquid crystal state with small anisotropy 
parameters (a = 1) will be quite small. In fact, the calcu- 
lation of energy differences between these states and the 
isotropic state are comparable to the estimated accuracy 
of our energy calculations. However, since the energy 
differences for larger a's show a definite tendency in all 
cases, we believe that the results are, significantly reliable 
(since the statistical uncertainty on any MC calculation 
is systematic, the energy differences may be even more 
accurate than the absolute energies). 

Tables I, II and III present the results for the calcula- 



tion of the LLL correlation energies obtained by means 
of Eqs. (11) or (13), using the angle averaged pair corre- 
lation functions (or static structure factors) for the three 
different forms of the interaction potential for a variety of 
quasi-2D layer widths A [see Eqs. (6,7)]. When A = all 
interaction potentials reduce to the Coulomb potential 
and in the case of the vzds {1^12) potential we note that 
b = 1.5/A. Results for filling factors 1^ = 1/3, 1/5 and 1/7 
of the LLL (for the potential Ui(ri2)) are also presented 
in Fig. 5. The results suggest that, in the LLL, for all 
the interaction potentials under consideration, a uniform 
liquid state is energetically more favorable than the liq- 
uid crystal state. For small values of a G (0,« 2], the 
liquid crystal states have an energy only slightly above 
the Laughlin liquid states (a = 0), however for larger a's 
this difference increases. 
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FIG. 5. Energy difference between ansiotropic states and the isotropic state {a = 0) AEa = Ea — Eq for filling factors 
f = 1/3, 1/5 and 1/7 in the LLL. These results correspond to the interaction potential vi{ri2) and are plotted as function of 
the quasi-2D layer thiclcness A. 
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FIG. 6. Energy difference between ansiotropic states and 
the isotropic state (a = 0) AEa = Ea — Eo for filling factor 
V = l/Z in the second excited valence LL (L = 2). These 
results correspond to the interaction potential form vi(ri2) 
and are plotted as function of the quasi-2D layer thickness A. 

Similar results are obtained in the first excited LL 
[L = 1 in Eq. (14), we omit the results for brevity]. For 
all forms of the interaction potential considered here, the 
correlation energy for anisotropic states is higher, once 



again leaving the Laughlin state stable. However, it is 

interesting to note that for the second excited LL {L = 2) 
the situation changes for the nematic states dX v = 1/3 
of the valence LL, where anisotropic states become en- 
ergetically favorable. Table IV shows the results for the 
energies, and energy differences, AE^^ = Ea — Eq 
(also shown in Fig. 6) between anisotropic states {a ^ 0) 
and the isotropic state (a = 0) for filling factor z/ = 1/3 
in the second excited LL {L = 2) obtained from poten- 
tial v\{r\2) (the results are quite similar for the other 
two forms of the potential). These results arc gener- 
ally consistent to what we found in the past using the 
hypernetted-chain (HNC) approximation [18,19]. 

A conclusion can be derived from the above results: 
generally speaking the isotropic states seem to be en- 
ergetically favorable, with the exception of the nematic 
state in the second excited LL. The explanation for this is 
simple: in the LLL the electron packets are simple gaus- 
sians, and it is clear that the best way to minimize their 
Coulomb repulsion is by placing the vortices responsible 
for the CF transformation [4, 5] precisely at the location 
of the electron themselves (a = 0). In higher LL's, the 
wavepackets take a more "ring-like" shape, and a finite 
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a permits a more optimal distribution of charge for the 
nematic case (but not for either the tetratic or hexatic). 



V. UNDERLYING CHARGE DENSITY WAVE IN 
THE ANISOTROPIC 2DOCP'S 

In view of the appearance of considerable density vari- 
ations in our MC simulations for larger values of a we in- 
vestigated the possible existence of an underlying CDW 
for the liquid crystalline states of Eq. (3). For this pur- 
pose it is useful to consider, once again, the 2D0CP 
analog system. Whereas considerable effort has been 
dedicated (and a consequent vast knowledge has been 
achieved) in the past to the treatment of the standard 
isotropic plasma (see e.g. Refs. [10,26,32,33]), little has 
been pursued for a system with anisotropic interactions, 
e.g. quadrupolar terms. 

Consider the classical distribution function (note: in 
this section we work in units of the magnetic length Iq): 



^i/ml^ oc e , where 
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where, as before, = ae*^^', 6^ = 27r(/i— l)/(m — 1), 
and /X € {1, 2, . . . , (m — 1)}. This potential energy corre- 
sponds to an "electrostatic potential" which is solution 
of a modified Poisson's equation 



m — 1 



V2[/3</,(r)] = -47r[p(r) + Pi^ - "/^)J + 47rmpo , (16) 



where = a(cos0^,sin0^). 
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and po = 1/ (27rm) is a neutralizing density. 



(17) 



nematic {u = 1/3) : ^ ^2^^ [1 + 2cos{aqx)] , 



Consider now the potential V generated by the addi- 
tion of some charge dp{r). This will cause a redistribution 
of the particles that form the plasma, inducing a density 
change [sec the discussion related to the definition of the 
pair correlation function, Eq. (8)]. 



Pindir) = J (fr' pob(r - r') - 1] 5p{v') , 



(18) 



The total charge, in reciprocal space, is therefore given 
by [see Eq. (9)]: 



Pt„t(k) = 5(k)(5p(k), 



leading to a total potential: 



am = 



471 5(k) 

A;2 



E 



(19) 



5p(k) . (20) 



This result neglects second order corrections in the dis- 
tribution functions and is, therefore, commonly referred 
to as the theory of linear screening. 

It is now interesting to investigate whether this po- 
tential allows for the formation of underlying CDW's in 
the 2D0CP. Assuming small variations from a uniform 
state, we allow for the particle density to vary from point 
to point according to: 



p{t) = Po + Pi cos(q • r) . 



(21) 



where q is the wavevector of the CDW and pi <C po- The 
0[pf] "excess energy" [10] per unit area is given by: 



Pi 



1 27r5'(q) 
2 



m— 1 



(22) 



It is evident that the charge neutrality sum rule 
(S'(q) oc for g ^ [10,25,26]) guarantees the ehm- 
ination of the singularity at g = leading to screening 
of the interaction. More interesting, however, is the fact 
that the excess energy becomes negative for a variety of 
wavevectors when a ^ 0. If we write Eq. (22) explicitly 
for the various states considered in this paper: 



tetratic {v = l/h) : ^ ^2*^ [1 + 2 cos{aqx) + 2 cos(agj^)] , 
2 q 

1 • / 1 /„N 1 27rS'(q).^ „ , s „ r / 1 - r / 1 'v/S MT 

hexatic [v = 1/7) : ^ — \^ + 2cos{aqx) + 2cos[a{--qx + —qy)\ + +2cos[a[--qx - —qy)\\ ■ 



(23) 



We can see that the most important configurations (those that make the potential PV minimum and max- 
imize their probability) correspond to charge density waves with wavevectors in the neighborhood of [34] 
aq/n ~ {(1, 0), (— 1. 0)} for the nematic, aq/n ~ {(1, 1), (1, — 1), (— 1, 1), (— 1, — 1)} for the tetratic, and aq/n ~ 
{(|, 0), (— |, 0), (|, ^), (|, — :^), (~|> (■"§' ~'^)} hexatic. This should produce a unidirectional CDW 

(a layered system, or smectic) underlying the nematic, with a characteristic wavelength A ~ 2a; a square lattice tilted 
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45° with lattice constant a ~ \/2a, and a triangular lattice with triangle side a = ^/Sa. Figure 7 depicts typical 
configurations during MC simulations with large a's. The characteristic CDW's have periods very close to those 
predicted above. 

One should note that these underlying CDW's arc extremely soft and fluctuations will render them invisible in the 
thermodynamic and ergodic limits. In our simulations, however, their effects are perceptible (see e.g. Fig. 1) for large 
values of the anisotropy parameter a because of phase locking at the boundaries. A detailed study of the fluctuations 
of these CDW's will be published elsewhere [35] . 
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FIG. 7. Typical electron configurations for a nematic {v = 1/3, a = 7, left panel), tetratic {v = 1/5, q = 8, center panel), 
and nematic {v = 1/7, a = 10, right panel). Note the formation of a CDW's with one, two and three different directors. 



VI. CONCLUSIONS 

In conclusion, wc have investigated the possibility of 
liquid crystal states in quasi-two-dimensional electron 
systems in strong magnetic flelds. We considered trans- 
lation invariant yet anisotropic states at filling factors 
1/ = 1/3, 1/5 and 1/7 of the lowest {L = 0), first ex- 
cited (L = 1) and second excited {L = 2) LL's. We 
found that the anisotropic states poscss an imdcrlying 
CDW along directors with the same symmetry group of 
the proposed state but these CDW's are "washed-out" 
by fiuctuations. Wc applied MC methods to calculate 
the (angle-dependent) pair correlation function and static 
structure factors for these states, which have permitted 
us to calculate the correlation energies for a variety of 
reasonable generalizations of the Coulomb potential that 
take into consideration the finite width of the quasi-2D 
layer. For all states and potentials under consideration 



the isotropic Laughlin state is found to be energetically 
favorable in the lowest and first excited LL, whereas we 
find an instability of the v = 1/3 nematic state in the 
second excited LL. 
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TABLE I. Correlation energy per particle in the LLL (in 
units of /do) for the liquid crystal (BRS) states at filling 

factor = 1/3 as a function of the anisotropy parameter a 
and quasi-2D layer width A. Three forms of the interaction 
potential were used. The three potentials reduce to the stan- 
dard Coulomb potential for A = 0. 



Interaction Ptrteutial: (■i(;-i2) 


a 


A = 0.0 


A = 0.5 


A = 1.0 


A = 1.5 


A = 2.0 


A = 2.5 


A = 3.0 





-0.4100 


-0.3362 


-0.2776 


-0.2327 


-0.1973 


-0.1700 


-0.1485 


1 


-0.4098 


-0.3353 


-0.2770 


-0.2319 


-0.1970 


-0.1698 


-0.1483 


2 


-0.3961 


-0.3234 


-0.2681 


-0.2257 


-0.1928 


-0.1669 


-0.1464 


3 


-0.3608 


-0.2926 


-0.2449 


-0.2093 


-0.1817 


-0.1597 


-0.1418 


4 


-0.3074 


-0.2435 


-0.2038 


-0.1763 


-0.1554 


-0.1387 


-0.1249 


Interaction Potential: V2(;rrz) 


a 


A = 0.0 


A = 0.5 


A = 1.0 


A = 1.5 


A = 2.0 


A = 2.5 


A = 3.0 





-0.4100 


-0.3286 


-0.2598 


-0.2107 


-0.1760 


-0.1507 


-0.1315 


1 


-0.4098 


-0.3277 


-0.2593 


-0.2104 


-0.1758 


-0.1505 


-0.1314 


2 


-0.3961 


-0.3162 


-0.2519 


-0.2058 


-0.1727 


-0.1483 


-0.1297 


3 


-0.3608 


-0.2859 


-0.2324 


-0.1936 


-0.1650 


-0.1433 


-0.1264 


4 


-0.3074 


-0.2370 


-0.1950 


-0.1661 


-0.1439 


-0.1265 


-0.1125 


Interaction Potential: vzds{^''i2) 


a 


A = 0.0 


A = 0.5 


A = 1.0 


A = 1.5 


A = 2.0 


A = 2.5 


A = 3.0 





-0.4100 


-0.3279 


-0.2748 


-0.2381 


-0.2112 


-0.1904 


-0.1738 


1 


-0.4098 


-0.3270 


-0.2741 


-0.2376 


-0.2107 


-0.1900 


-0.1735 


2 


-0.3961 


-0.3160 


-0.2657 


-0.2310 


-0.2053 


-0.1854 


-0.1696 


3 


-0.3608 


-0.2873 


-0.2439 


-0.2138 


-0.1914 


-0.1739 


-0.1597 


4 


-0.3074 


-0.2408 


-0.2054 


-0.1813 


-0.1633 


-0.1491 


-0.1375 



TABLE II. Correlation energy per particle in the LLL (in 
units of /do) for the liquid crystal (BRS) states at filling 

factor = 1/5 as a function of the anisotropy parameter a 
and quasi-2D layer width A. Three forms of the interaction 
potential were used. The three potentials reduce to the stan- 
dard Coulomb potential for A = 0. 



Interaction Ptrteutial: ('i(;-i2) 


a 


A = 0.0 


A = 0.5 


A = 1.0 


A = 1.5 


A = 2.0 


A = 2.5 


A = 3.0 





-0.3274 


-0.2811 


-0.2420 


-0.2094 


-0.1825 


-0.1603 


-0.1419 


1 


-0.3273 


-0.2810 


-0.2419 


-0.2094 


-0.1825 


-0.1603 


-0.1419 


2 


-0.3265 


-0.2803 


-0.2413 


-0.2089 


-0.1821 


-0.1600 


-0.1418 


3 


-0.3121 


-0.2674 


-0.2312 


-0.2014 


-0.1767 


-0.1563 


-0.1392 


4 


-0.2775 


-0.2362 


-0.2064 


-0.1829 


-0.1635 


-0.1472 


-0.1333 


5 


-0.2216 


-0.1836 


-0.1601 


-0.1432 


-0.1296 


-0.1181 


-0.1081 


Interaction Potential: V2{r\_2) 


a 


A = 0.0 


A = 0.5 


A = 1.0 


A = 1.5 


A = 2.0 


A = 2.5 


A = 3.0 





-0.3274 


-0.2743 


-0.2365 


-0.2086 


-0.1873 


-0.1704 


-0.1566 


1 


-0.3273 


-0.2743 


-0.2365 


-0.2086 


-0.1873 


-0.1704 


-0.1566 


2 


-0.3265 


-0.2767 


-0.2303 


-0.1928 


-0.1641 


-0.1422 


-0.1251 


3 


-0.3121 


-0.2639 


-0.2215 


-0.1870 


-0.1603 


-0.1396 


-0.1233 


4 


-0.2775 


-0.2329 


-0.1997 


-0.1730 


-0.1513 


-0.1338 


-0.1196 


5 


-0.2216 


-0.1801 


-0.1561 


-0.1378 


-0.1222 


-0.1088 


-0.0975 


Interaction Potential: vzDs{ri2) 


a 


A = 0.0 


A = 0.5 


A = 1.0 


A = 1.5 


A = 2.0 


A = 2.5 


A = 3.0 





-0.3274 


-0.2743 


-0.2365 


-0.2086 


-0.1873 


-0.1704 


-0.1566 


1 


-0.3273 


-0.2743 


-0.2365 


-0.2086 


-0.1873 


-0.1704 


-0.1566 


2 


-0.3265 


-0.2736 


-0.2359 


-0.2082 


-0.1869 


-0.1701 


-0.1563 


3 


-0.3121 


-0.2615 


-0.2265 


-0.2006 


-0.1807 


-0.1648 


-0.1519 


4 


-0.2775 


-0.2324 


-0.2037 


-0.1825 


-0.1659 


-0.1525 


-0.1414 


5 


-0.2216 


-0.1819 


-0.1598 


-0.1439 


-0.1313 


-0.1211 


-0.1124 
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TABLE III. Correlation energy per particle in the LLL (in 
units of /do) for the liquid crystal (BRS) states at filling 

factor = 1/7 as a function of the anisotropy parameter a 
and quasi-2D layer width A. Three forms of the interaction 
potential were used. The three potentials reduce to the stan- 
dard Coulomb potential for A = 0. 



Interaction Ptrteutial: (■i(;-i2) 


a 


A = 0.0 


A = 0.5 


A = 1.0 


A = 1.5 


A = 2.0 


A = 2.5 


A = 3.0 





-0.2827 


-0.2491 


-0.2198 


-0.1944 


-0.1727 


-0.1541 


-0.1383 


1 


-0.2827 


-0.2491 


-0.2198 


-0.1944 


-0.1727 


-0.1541 


-0.1383 


2 


-0.2826 


-0.2491 


-0.2198 


-0.1944 


-0.1727 


-0.1541 


-0.1383 


3 


-0.2807 


-0.2473 


-0.2184 


-0.1933 


-0.1719 


-0.1536 


-0.1379 


4 


-0.2492 


-0.2185 


-0.1945 


-0.1745 


-0.1573 


-0.1425 


-0.1296 


5 


-0.1917 


-0.1643 


-0.1467 


-0.1334 


-0.1223 


-0.1124 


-0.1035 


Interaction Potential: V2{r\_2) 


a 


A = 0.0 


A = 0.5 


A = 1.0 


A = 1.5 


A = 2.0 


A = 2.5 


A = 3.0 





-0.2827 


-0.2470 


-0.2123 


-0.1821 


-0.1576 


-0.1382 


-0.1227 


1 


-0.2827 


-0.2470 


-0.2124 


-0.1822 


-0.1577 


-0.1383 


-0.1227 


2 


-0.2826 


-0.2470 


-0.2123 


-0.1821 


-0.1576 


-0.1382 


-0.1227 


3 


-0.2807 


-0.2452 


-0.2110 


-0.1813 


-0.1571 


-0.1378 


-0.1225 


4 


-0.2492 


-0.2164 


-0.1894 


-0.1659 


-0.1460 


-0.1297 


-0.1162 


5 


-0.1917 


-0.1621 


-0.1443 


-0.1295 


-0.1160 


-0.1040 


-0.0936 


Interaction Potential: vzDsiru) 


a 


A = 0.0 


A = 0.5 


A = 1.0 


A = 1.5 


A = 2.0 


A = 2.5 


A = 3.0 





-0.2827 


-0.2436 


-0.2141 


-0.1914 


-0.1735 


-0.1591 


-0.1471 


1 


-0.2827 


-0.2437 


-0.2141 


-0.1915 


-0.1736 


-0.1591 


-0.1472 


2 


-0.2826 


-0.2436 


-0.2141 


-0.1914 


-0.1735 


-0.1591 


-0.1471 


3 


-0.2807 


-0.2419 


-0.2128 


-0.1904 


-0.1727 


-0.1584 


-0.1465 


4 


-0.2492 


-0.2147 


-0.1906 


-0.1720 


-0.1572 


-0.1450 


-0.1348 


5 


-0.1917 


-0.1627 


-0.1454 


-0.1323 


-0.1216 


-0.1126 


-0.1050 



TABLE IV. Correlation energy per particle in the second 
excited LL, L = 2, (in units of e^/eZo) for the liquid crys- 
tal (BRS) states at filling factor = 1/3 as a function of 
the anisotropy parameter a and quasi-2D layer width A. The 
form v\{r\2) for the interaction potential was used. 



Interaction Potential: «i(ri2) 



a 


A = 0.0 


A = 0.5 


A = 1.0 


A = 1.5 


A = 2.0 


A = 2.5 


A = 3.0 





-0.2642 


-0.2139 


-0.1872 


-0.1662 


-0.1485 


-0.1335 


-0.1207 


1 


-0.2653 


-0.2146 


-0.1875 


-0.1663 


-0.1486 


-0.1335 


-0.1208 


2 


-0.2693 


-0.2169 


-0.1881 


-0.1663 


-0.1483 


-0.1333 


-0.1206 


3 


-0.2708 


-0.2158 


-0.1852 


-0.1631 


-0.1455 


-0.1310 


-0.1188 
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